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The second order correction to free energy due to the interaction between elec- 
trons is calculated for a quasi-one-dimensional conductor exposed to a magnetic field 
perpendicular to the chains. It is found that specific heat, magnetization and torque 
oscillate when the magnetic field is rotated in the plane perpendicular to the chains 
or when the magnitude of magnetic filed is changed. This new mechanism of ther- 
modynamic magnetic oscillations in metals, which is not related to the presence of 
any closed electron orbits, is applied to explain behavior of the organic conductor 
(TMTSF) 2 C10 4 . 
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In a recent study [|TJ of the quasi-one-dimensional conductor (TMTSF) 2 C104, a new 
physical effect has been found. In this experiment, a magnetic field was rotated in the plane 
perpendicular to the conducting chains and the torque along the chains was measured. The 
torque, as a function of the angle 9 between the magnetic field and the c-axis, was observed 
to exhibit oscillations approximately periodic in tan#. Since the torque is a derivative of the 
free energy F with respect to 9, then F must also oscillate upon rotation of the magnetic 
field. In earlier experiments, magnetization and specific heat || were measured to be 
oscillating functions of the magnetic field magnitude (the so called "fast" oscillations). Such 
angular and fast thermodynamic magnetic oscillations were totally unexpected theoretically 
for a quasi-one-dimensional metal. In momentum space, an electron in a magnetic field 
moves along an orbit obtained as the intersection of a constant energy surface and a plane 
perpendicular to the magnetic field. Magnetic oscillations of the free energy (de Haas- 
van Alfven effect) are known in the situation when the orbits at the Fermi surface are 
closed. In this case, the orbits are quantized, the energy spectrum is discrete and the 
oscillations in F are due to the Landau levels crossing the Fermi level. However, in a 
quasi-one-dimensional conductor, the Fermi surface consists of two slightly warped planes 
perpendicular to the direction of the chains. Thus, if a magnetic field has been applied 
perpendicular to the chains, all electron orbits at the Fermi surface are open, the energy 
spectrum is continuous and no oscillations in F are expected for non-interacting electrons 
||]. In the present paper, the lowest order correction to F due to the interaction between 
electrons is calculated for a quasi-one- dimensional conductor and is shown to exhibit angular 
and fast magnetic oscillations. 

The phenomenon of angular oscillations in a quasi-one- dimensional conductor was pre- 
dicted theoretically by Lebed |J who calculated the oscillations of the magnetic-field-induced 
spin-density-wave (FISDW) transition temperature. Experiment || did not reveal these os- 
cillations, however, the studied range of angles might not have been big enough to see the first 
oscillation [jT]. The angular oscillations in the conductivity were predicted ]7J and observed 
f|J^Jl|]- Two other models for the oscillations in the conductivity were later suggested |10 



None of mentioned theoretical papers considered the oscillations of free energy. Recently, 
the contribution of the FISDW order parameter fluctuations to F has been calculated and 



been shown to exhibit angular oscillations [11]. In the present paper, the approach of Ref. 



H and J7J is used to calculate the lowest order correction to F. This theory is applicable in 
the metallic phase sufficiently far from the phase boundary, while theory [|IT[] is applicable 
in the vicinity of the phase transition to the FISDW state. 

Let us consider a square lattice of parallel chains with a spacing b between the chains, 
x-axis parallel the chains. The electron hopping integrals between the nearest chains in the 
y and z directions are equal to tb and t c , respectively; t c <C tb <C Ep, where Ep is the Fermi 
energy. Magnetic field H is applied in the (y, z) plane at an angle 9 measured from z-axis. 



The electron Hamiltonian without the interaction term can be written as 12 



H aa = —aivpdx + 2tbCos(k y — ebHx cos 6*/c) 

+2t c cos(k z — ebHx sin 9/ c) — fisHa, (1) 

where the index a = ± labels electrons whose momenta along the chains are close to ±/c_f, 
vf and kp are the Fermi velocity and the Fermi momentum, k y and k z are momenta along 
the y- and z-axes, a is the spin index, /i# is the Bohr magneton, e is the electron charge, 
c is the velocity of light and h = 1. The two cosine terms in (|l|), which describe electron 
motion perpendicular to the chains, become in magnetic field two periodic potentials along 
the chains. Their periods are commensurate at rational values of tan#. This effect is 
responsible for the angular oscillations of free energy. 

The lowest order non-trivial correction to the free energy AF due to the interaction 
between electrons is given by the diagram shown in Fig. 1. Using Green functions, derived 
from Eq. ([!]) [12|], one can find the following expression for the free energy per electron: 



27rVT 3 f°° , cosh(rx) „ . . „ , , 

A^ = -r^r dx , ) \ f y (x)f z (x), (2) 

ep^l Jo smh (tx) 

f y (x) = (jQ(Vf,sin(xcos6 l ) sinfc/ cos9))k, (3) 
fz{%) — (</o(Ksin(a;sin#) sin /c/ sin 6*))^. (4) 



Here Ef = vf^f, T is the temperature, = ebvFH/cis the characteristic energy of magnetic 
field, t = A-kT/VL, Vu c \ = St b ^/Q, J is the Bessel function and (...)k means the averaging 
over k. g 2 = {g2~ 9i) 2 + where g\ and #2 are dimensionless interaction constants between 
electrons, normalized in such a way that the mean-field SDW transition temperature is 
proportional to exp(— l/^)- Another contribution to AF, which comes from spin reversal 
scattering, contains an additional factor cos(4fiBHx/Q) under the integral in Eq. (§) and 
g\ appears instead of g 2 . This contribution is more complex and will not be studied here. 

Integral @) is divergent at small x. However, observable quantities considered below are 
given by convergent integrals. The correction to the specific heat per electron, derived from 
Eq. d2|), takes the form: 

Q /*00 / \ ' r~ r~ 

where the primes denote the third derivative with respect to x and C = tt 2 T/3ef is the 
specific heat of non-interacting electrons. In Fig. 2, AC (8) is shown for three different 
values of r (0.1, 0.2 and 0.4) at H = 20 and V c = 2. AC exhibits sharp peaks when tan.9 
attains rational values, which is where the periods of the functions f y {x) and f z {x) in ([S]) 
become commensurate. The widths of the peaks go to zero when T — > and the curve 
becomes fractal. Each peak is split because the term [x 2 j sinh 2 x)'" in (H) oscillates once 
over the range [0, oo]. This splitting vanishes as T — > 0. It follows from (||) and (f|) that the 
transformation 9 — > 45° — 9 interchanges Vj, and V c in (|5|), having small effect on the plots 
AC (9) in Fig. 2, despite the employed values of V& and V c differ by an order of magnitude. 
This means that AC does not depend qualitatively on the precise values of Vb and V c if they 
are greater than 2. According to Ref. JlJ, in (TMTSF) 2 C10 4 , Q./H = 1.8 K/T, but due to 
the unit cell doubling the two times bigger value should be used. Thus, the field e.g. 5 T 
corresponds to Q = 18 K and the values of the parameters r and V in Fig. 2 correspond to 
T = 0.14K, 0.29K, 0.57K, t b = 45K and t c = 4.5K. The dependence of specific heat on the 
orientation of the magnetic field has not been measured yet experimentally. 

The correction to the magnetization per electron along the z-axis AM Z = —dAF/dH z , 
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where H z = HcosO, is given by the following expression: 

qHh^l f°° coshf-^Xxcosx — sinx) 
AM, = — — / dx- Kcos0jK ' 



2e F Hh ( Sinh (^)^)3 

, T M sin x sin k . T ,V b sin x sin fc x . 
sm/cJi J )) k f z { -), 6 

cos 9 COS 9 COS 9 

Unlike AC, AM z (0) has a limit of T = which is shown in Fig. 3 at V h = 20 and V c = 2. Note 
that the magnetic energy H Z AM Z may be much greater than the Zeeman energy (^bH) 2 /ef 
and the orbital energy (t 2 cos 2 9 + t 2 sin 2 9)Q 2 /12e 3 F of the non-interacting electrons (the 
latter formula is a result beyond quasi-classics) . The magnetization is found to be negative 
(diamagnetic) in agreement with the experiments in the high field reentrant metallic phase 
of (TMTSF) 2 C10 4 0. The torque is determined by the expression N x = — AM Z H cos9 + 
AM y H sin9, where AM y is given by Eq. (||) after interchanging tb with t c and cos 9 with 
sin#. The torque is plotted in Fig. 3 at r = and r = 0.25 together with an experimental 



curve taken at T = 0.4 K, H = 5 T ||14j| . Using the cited above value of Q/H, one finds 
t = 0.28 in the experiment, close to the value employed for the theoretical curve c. There 
is an overall qualitative resemblance between the theoretical and the experimental curves. 
For a detailed comparison one should take into account the different spacings of the chains 
along the y- and z-axes and the triclinic structure of the real material. Note that, with few 
exceptions, the peaks in Fig. 3 appear when tan# is integer, not fractional as in Fig. 2, 
in agreement with the experimental observation This happens because at r = the 
aperiodic function of x under the integral decreases rapidly when x > 2tt, while the last 
function in (^) has the (longer) period of the order of ir/ tan#, thus the resonances cannot 
take place at small tan#. 

In Fig. 4, AM Z (||) at 9 = is plotted vs the variable l/H? which is proportional to 
the magnetic field. AM Z exhibits oscillations periodic in inverse magnetic field, which may 
correspond to the fast oscillations of magnetization observed in Ref . 0] . The period of the 
oscillations A(l/H) is determined by the condition AVb = it and equals irebvF/8ctb. The 
numbers of oscillations in Fig. 4 are from 3 to 15. Using the experimental value 255T for 
the frequency of the fast oscillations || one finds tb = 180K. A similar mechanism for the 



fast oscillations of conductivity was studied in Ref. | 15| . However, we have not found these 
oscillations in specific heat (W). 



Another possible explanation of the fast oscillations ||T6[ takes into account that the 
unit cell of (TMTSF) 2 C104 is doubled in the y direction due to a crystal superstructure, 
thus the energy of a chain staggers by the value ±k along the y-axis. Suppose, that k ^> 
tb- In this case, the conduction band is split into two bands with the Fermi momenta 
kp ± k/vf, and the effective transverse hopping integral tl is equal to t\j2K. In one band, 
the electron wave functions are predominant on even chains; in another band, on odd chains. 
Interaction between the electrons belonging to the same band reproduces the formulas given 
above with t\ substituted for % and the spacing along y-axis doubled. This prescription 
immediately explains why in experiments [|l]|| the angular oscillations were found only at 
even, instead of all, integer values of tan 9 defined with respect to the original lattice without 
superstructure. In addition, the interaction of the electrons belonging to the different bands 
gives the following contribution to the free energy (at 8 = 0): 

2TT 3 g 2 T 3 tt , cosh(xr') . t2 . ta . .. 

A F = TTTT^ dx ■ J, M (V c xsmp)) p 

EpiiK, Jo smh (XT') 

([(1 + 2 cos 2 x) Jq (V b sin x sin k) — Jf(V b * sinxsinfc)] 
sin 2 k cos(2_Rx) + 4J Q (V b * sin x sin k)Ji(V b * sin x sin k) 

sin /c cosx sin(2i?x))fc, (7) 

where V b * = 4t b /Q, t' = r/2 and R = 2k/Q. The correction to the specific heat A'C is 
obtained from (0) in the same way as Eq. ([5D is obtained from Eq. (0) and is plotted vs R 
in Fig. 5. A'C exhibits oscillations periodic in inverse magnetic field which may correspond 
to the fast oscillations of specific heat found in Ref. 0. The period of the oscillations is 
determined by the condition AR = 1 and is given by A(l/H) = ebvp/2cK, which corrects 
by a factor of 2 the expression found in Ref. [[16]. The numbers of oscillations in Fig. 5 are 
from 2 to 40. Using the experimental frequency of the fast oscillations 257T |3[] one finds 
2k = 460K. Analogous oscillations should also exist in magnetization, which can be derived 
from Eq. (|7|) as well. If, on the other hand, the opposite limit k, < t b holds, then, in order 



to explain why the odd peaks are missing in experiments it is necessary to conclude 
that the field H z at H = 5 T is below the magnetic breakdown field H = K^c/t^ebvp. In 
this case, the mechanism described above still works qualitatively but with the frequency of 
the fast oscillations proportional to £& instead of k. If the odd peaks appear experimentally 
in higher fields, this can give the estimates of H and k. 

Fast oscillations were also observed in the material (TMTSF^PFg, which does not have 
an anion superstructure fl7j • This fact supports the first explanation of the fast oscillations 
which does not invoke the superstructure. On the other hand, these oscillations were ob- 
served thus far only in the FISDW phase in transport measurements. Thus, formally, the 
fast oscillations in (TMTSF) 2 PF 6 are outside the applicability range of presented theory 
and a definitive choice between two explanations cannot be made (see also Ref. [fL8[|). 

In conclusion, the second order correction to the free energy due to the interaction 
between electrons has been calculated for a quasi-one- dimensional metal. Angular and fast 
magnetic oscillations of specific heat, magnetization and torque were found in qualitative 
agreement with experiment. Two different mechanisms of the fast oscillations, related and 
unrelated to the specific crystal structure of (TMTSF) 2 C10 4 , were observed. Unlike in the 
standard theory of metals, the oscillations are completely due to the interactions between 
electrons and cannot be interpreted in terms of the presence of some closed electron orbits. 
The oscillations involve the energy scales of the cyclotron frequency Q, the hopping integral 
tf) or the anion superstructure splitting k, which are much higher than the characteristic 
energy scale of the deviation from nesting and the FISDW transition temperature. That 
explains why the oscillations persist in the FISDW phase as observed experimentally. 

The author thanks M. J. Naughton for sending torque data shown in Fig. 3; J. S. Brooks, 
who suggested to calculate specific heat; both of them and P. M. Chaikin, for discussions. 
This work was supported by the NSF Grant No. DMR 89-06958. 



7 



REFERENCES 

* On leave from L. D. Landau Institute for Theoretical Physics, Moscow, 117940, Russia. 

1 M. J. Naughton et al, Phys. Rev. Lett. 67, 3712 (1991). 

2 X. Yan et al, Phys. Rev. B 36, 1799 (1987); R. V. Chamberlin et al., Phys. Rev. Lett. 
60, 1189 (1988). 

3 N. A. Fortune et al, Phys. Rev. Lett. 64, 2054 (1990). 

4 I. M. Lifshits, M. Ya. Azbel and M. I. Kaganov, Electron Theory of Metals (Consultants 
Bureau, NY., 1973). 

5 A. G. Lebed, Pis'ma Zh. Exp. Teor. Fiz. 43, 137 (1986) [JETP Lett. 43, 174 (1986)]. 

6 X. Yan et al., Solid State Commun. 66, 905 (1988). 

7 A. Lebed and P. Bak, Phys. Rev. Lett. 63, 1315 (1989). 

8 G. S. Boebinger et al., Phys. Rev. Lett. 64, 591 (1990); M. J. Naughton et al, Mater. 
Res. Soc. Symp. Proc. 173, 257 (1990). 

9 T. Osada et al., Phys. Rev. Lett. 66, 1525 (1991). 

T. Osada, S. Kagoshima and N. Miura, Technical report of ISSP, Univ. of Tokyo, Ser. A, 
No. 2470 (1991); K. Maki, Phys. Rev. B 45, 5111 (1992). 

1 A. Bjelic and K. Maki, Phys. Rev. B, to be published. 

2 L. P. Gor'kov and A. G. Lebed', J. Phys. (Paris) Lett. 45, L433 (1984). 

3 V. M. Yakovenko, Zh. Exp. Teor. Fiz. 93, 627 (1987) [Sov. Phys. JETP 66, 355 (1987)]. 

4 The author is grateful to M. J. Naughton for permission to publish the raw experimental 
torque data of Ref. [[[]]. The curve d in Fig. 3 actually represents the data between two 
main minima in the angle interval from 6° to —90° mapped for the purpose of qualitative 
comparison to the interval 0° to 90°. 

8 



K. Yamaji, J. Phys. Soc. Jpn. 55, 1424 (1986). 

S. A. Brazovskii and V. M. Yakovenko, Pis'ma Zh. Exp. Teor. Fiz. 43, 102 (1986) [JETP 
Lett. 43, 134 (1986)]. 

J. P. Ulmet et al, J. Phys. (Paris) Lett. 46, L535 (1985); J. R. Cooper et al, Phys. Rev. 
Lett. 63, 1984 (1988); S. T. Hannahs et al., ibid. 63, 1988 (1988); W. Kang et al., Phys. 
Rev. B, to be published. 

A. G. Lebed', Phys. Scripta 39, 386 (1991). 



9 



FIGURES 

FIG. 1. The second order correction to the free energy. The lines, labeled by + and -, represent 
Green functions of the electrons with momenta close to ±&p. The vertices represent the amplitudes 
of interaction between electrons. 

FIG. 2. Normalized correction to specific heat —AC /3g 2 Co vs magnetic field rotation angle 
8. Parameter r is equal to 0.1 (curve a), 0.2 (b) and 0.4 (c). The arrows at the top indicate the 
rational values of tan8. 

FIG. 3. Normalized magnetization — 16AM z Hsf /g 2 Q 2 (curve a) and normalized torque 
lQN x £p/g 2 Q 2 (curve b, t = 0, and curve c, r = 0.25) calculated as functions of 9. The curve 
d represents the experimental torque data due to M. J. Naughton 14]. 



FIG. 4. Normalized magnetization — 1QAM z Hef / g 2 Q 2 at T = vs normalized magnetic field 
l/Vb = ebvpH/8tbc; t c = in curve a and t c = O.ltb in curve b. 

FIG. 5. Additional normalized correction to specific heat in the presence of the special crystal 
superstructure —A'CK 4 /6g 2 t^Co vs normalized inverse magnetic field R = 2kc/ ebvpH . Parameters 
are: t* b = k/8, t c = k/80 and T = k/20tt 2 . 
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